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ABSTRACT — Allozymic variation of Japanese dace Tribolodon hakonensis was examined at 20 loci in 
17 populations derived from Hokkaido to Kyushu. Intrapopulational variability indicated by the 
proportion of loci polymorphic and heterozygosity was close to the standard levels in vertebrates. On 
the other hand, allelic frequencies at five loci of the twenty were remarkably different among 
populations, and even replacement of predominant alleles was observed. The extent of interpopula- 
tional differentiation was estimated by coefficient of gene differentiation (G ST ) and genetic distance 
(D). Levels of these indices were remarkably high compared with the levels known among popula- 
tions in many vertebrates. Based on a dendrogram, we classified the populations into three local 
groups; 1) Northern Group, 2) Southern Group, and 3) Lake Biwa Group. These results suggest that 
the populations of different local groups have been highly differentiated from each other and the 
differentiation is due to geographical isolation. 


INTRODUCTION 

The genus Tribolodon which belongs to the 
family Cyprinidae has a unique characteristic on 
distribution: it is widely distributed from the upper 
reaches of rivers to the coastal regions of the sea in 
Japan and adjacent countries. The other genera of 
Cyprinidae do not show such a wide distribution. 
Therefore fishes of Tribolodon have been exten- 
sively studied from the viewpoints of physiology 
[1] and ecology [2]. However, their taxonomy had 
been confused and greatly different among tax- 
onomists [3-5], because the species of this genus 
have not been morphologically differentiated 
enough from each other. Our previous study using 
allozyme markers clearly showed that the four 
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species of Tribolodon classified by Nakamura [4] 
had been genetically differentiated enough as 
highly as the interspecific levels in vertebrates [6, 
7]. Thus, biochemical analysis provided good 
information as for taxonomy. In addition, we 
examined T. hakonensis from the waters in 
Fukushima Pref. and Kochi Pref. in the previous 
work, and found the possibility that this species 
contained the highly differentiated populations [7]. 

In this work, we further analysed T. hakonensis 
collected at more locations from Hokkaido to 
Kyushu, evaluated the level of intraspecific 
variability and the extent of interpopulational 
differentiation, and discussed the evolutionary 
aspects in populations of T. hakonensis. 

MATERIALS AND METHODS 

Specimens 

We collected a total of 883 specimens of T. 
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hakonensis in 17 waters from Hokkaido to 
Kyushu. Sample locations are shown in Figure 1, 
and sample sizes are given in Table 3. Fish were 
caught at a single location within a few days by a 
gill net, a casting net, and angling, and were frozen 
by dry ice at once and stored under — 20°C. In the 
Otsuchi, the Monobe, and the Shimanto Rivers, 
specimens were collected at a few locations, the 
upper reaches, and the lower reaches or the coastal 
regions of the sea, which were isolated by dums at 
present. 

Electrophoretic analyses 

Horizontal starch gel electrophoresis was per- 
formed upon drip from tissue samples as described 
by Numachi [8]. The following four buffer systems 
(Table 1) were used: CAPM (citric acid, 4-(3- 
aminopropyl)-morpholine), pH 6.0 described by 
Clayton and Tretiak [9]; CAEA (citric acid, 


N-(3-aminopropyl)-diethanol-amine), pH 7.0 and 
CT (citric acid, tris), pH 8.0 described by Numachi 
et al. [10]; and TBE (tris, boric acid, EDTA), pH 
8.7 described by Numachi [11]. The staining 
methods were according to Shaw and Prasad [12], 
and Taniguchi and Numachi [13]. The gels stained 
were washed and dried by the method of Numachi 
[14]. A list of proteins analysed, their abbrevia- 
tions, E. C. numbers, locus designations, tissues, 
and buffers employed are shown in Table 1. 
Tissues of skeletal muscle, liver, heart, and blood 
were used. Mitochondrial and supernatant iso- 
zymes of AAT, IDH, MDH, and ME were 
discriminated. 

Locus and allele designation followed Hanzawa 
et al. [15]. The most common alleles in the 
populations of the Northern Group were desig- 
nated as 100 or —100, and the other alleles were 
designated by the relative differences in elec- 



Fig. 1. The sites where specimens of T. hakonensis were sampled. 
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Table 1. Enzymatic and non-enzymatic proteins 
electrophoretic buffers employed 

examined, locus 

designations. 

tissues, and 

Protein (Abbreviation; E.C. Number) 

Locus 

Tissue 

Buffer 

Aspartate aminotransferase (A AT; 2.6. 1.1) 

m-Aat 

muscle 

CT 

Alcohol dehydrogenase (ADH; 1.1. 1.1) 

Adh 

liver 

CAEA 

^-Glycerophosphate dehydrogenase (a-GDH; 1.1. 1.8) 

a-Gdh-2 

liver 

CAEA 

Glucosephosphate isomerase (GPI; 5.3. 1.9) 

Gpi-1 

heart 

TBE 


Gpi-2 

heart 

TBE 

Hemoglobin (HB) 

Hb-1 

blood 

CT 


Hh-2 

blood 

CT 

Isocitrate dehydrogenase (IDH; 1.1.1.42) 

s-Icd-1 

liver 

CAEA 


s-lcd-2 

liver 

CAEA 

Lactate dehydrogenase (LDH; 1.1.1.27) 

Ldh-1 

heart 

CT 


Ldh-2 

muscle 

CT 


Ldh-3 

liver 

CT 

Malate dehydrogenase (MDH; 1.1.1:37) 

s-Mdh-1 

muscle 

CAPM 


s-Mdh-2 

liver 

CAPM 


m-Mdh 

muscle 

CAEA 

Malic enzyme (ME; 1.1.1.40) 

s-Me 

muscle 

CAEA 

Phosphoglucomutase (PGM; 2. 7. 5.1) 

Pgm 

muscle 

CAEA 

6-Phosphogluconate dehydrogenase (6-PGD; 1.1.1.44) 

6-Pgd 

liver 

CAEA 

Sarcoplasmic protein (SP) 

Sp-2 

muscle 

CAEA 


Sp-3 

muscle 

CAEA 


trophoretic mobility of respective gene products 
with the reference allele, 100 or —100. 

Statistical analyses 

Amount of intrapopulational variability was 
evaluated by the proportion of loci polymorphic 
and heterozygosity. Degree of interpopulational 
differentiation was evaluated by coefficient of gene 
differentiation (G sr ) [16] and genetic distance (D) 
[17]. G st value, which is calculated based on 
average heterozygosity, indicates the degree of 
allelic differentiation over all populations. D 
value, which is calculated based on allelic frequen- 
cies, indicates the degree of differentiation be- 
tween pairs of populations. Divergence time 
between the populations was estimated based on D 
value and gene substitution rate of 10 -7 [16]. 

RESULTS 

Intrapopulational variation 

Twenty loci controlling 12 kinds of enzymatic 


and nonenzymatic proteins were examined. The 
genetic control by 20 loci was mainly deduced from 
the phenotypes of different species and their 
hybrids [6, 7, 15]. Allozymic variation was found 
at 14 of the 20 loci in sample populations of T. 
hakonensis. Particularly, various polymorphisms 
were observed at 5 loci, Adh , a-Gdh-2 , Gpi-2 , 
s-Mdh-2 , 6-Pgd, and each allelic control was 
interpreted (Fig. 2). Deviations of observed num- 
ber of phenotypes from their expectations under 
Hardy-Weinberg equilibrium were not significant 
in all sample populations except the upper reach of 
the Shimanto River (at 6-Pgd, P<0.05). 

Intrapopulational variability of T. hakonensis 
indicated by proportion of loci polymorphic and 
heterozygosity was given in Table 2. The propor- 
tion of loci polymorphic ranged from 0.05 to 0.50, 
and the mean value indicated 0.22 + 0.10. Heter- 
ozygosity observed ranged from 0.019 to 0.085, 
and the mean value indicated 0.041+0.019. 
Levels of these two indices were similar among 
sample populations. Heterozygosity observed 
agrees with its expectation (H°/H E =1). 
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Table 2. Estimates of genetic variability at 20 loci within populations of T. hakonensis 


Location 

Proportion of 
loci polymorphic 

Heterozygosity 

h°/h e 

(Observed; H°) 

(Expected; H E ) 

Abira R. 

0.30 

0.039 

0.040 

0.976 

Otsuchi 





upper reach 

0.20 

0.019 

0.020 

0.950 

Otsuchi Bay 

0.30 

0.028 

0.028 

1.000 

Hirose R. 

0.25 

0.046 

0.047 

0.979 

Kinu R. 

0.50 

0.069 

0.075 

0.922 

Kano R. 

0.30 

0.072 

0.081 

0.889 

Aga R. 

0.25 

0.035 

0.037 

0.952 

Iburibashi R. 

0.15 

0.030 

0.029 

1.039 

Yoshii R. 

0.10 

0.019 

0.019 

0.947 

Sufu R. 

0.20 

0.023 

0.021 

1.098 

Chikugo R. 

0.15 

0.033 

0.033 

1.014 

L. Biwa 

0.40 

0.085 

0.089 

0.954 

Kumano R. 

0.25 

0.072 

0.074 

0.970 

None R. 
Monobe R. 

0.15 

0.034 

0.030 

1.125 

upper reach 

0.05 

0.026 

0.024 

1.111 

lower reach 

0.15 

0.029 

0.024 

1.212 

Kure R. 

0.10 

0.028 

0.031 

0.917 

Shimanto R. 





upper reach 

0.25 

0.050 

0.059 

0.848 

middle reach 

0.20 

0.048 

0.049 

0.974 

lower reach 

0.20 

0.042 

0.042 

1.000 

Hitotsuse R. 

0.20 

0.034 

0.034 

1 .003 


(0.22 + 0.10) 

(0.041 + 0.019) 




Inter pop ulational differentiation 

Gene constitution of sample populations col- 
lected at 21 locations in 17 waters from Hokkaido 
to Kyushu were compared (Table 3). Allelic 
frequencies among subpopulations within the same 
waters were very similar, and a significant devia- 
tion among allelic frequencies was not observed at 
all loci except Adh and a-Gdh-2 in the Shimanto 
River population. These subpopulations were 
isolated by dums at present, and fish of the Otsuchi 
Bay was particularly the amphidromous type. 
However, gene constitution of these subpopula- 
tions within the same waters was thus very similar. 

On the other hand, remarkable differences 
among populations derived from the different 
waters were found at five loci, Adh,a-Gdh-2, 


Gpi-2 , s-Mdh-2, 6-Pgd. At Adh and a-Gdh-2 of 
the five, predominant alleles of 100 in the popula- 
tions of the Pacific coast of the northern Japan and 
of the coast of the Sea of Japan were replaced with 
the other alleles (119 or 0) in those of the Pacific 
coast of the southern Japan. Similarly, predomi- 
nant allele of 100 at 6-Pgd was replaced with 127 
only in the Lake Biwa and the Kumano River 
populations. Population specific alleles with high 
frequencies were found, such as 176 at Adh in the 
Lake Biwa population and —153 at Gpi-2 in the 
Kumano River population. 

The extent of interpopulational differentiation 
was evaluated by two kinds of indices, coefficient 
of gene differentiation (G ST ) and genetic distance 
(D). As for the Otsuchi, the Monobe, and the 
Shimanto River systems, specimens collected at 


Table 3. Allelic frequencies for 14 loci in populations of T. hakonensis 
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Table 4. Estimates of Nei’s genetic distance between pairs of 17 populations of T. hakonensis based on 
allelic frequencies at 20 loci 




17 

16 

15 

14 

13 12 11 

10 

9 

8 

7 

6 

5 

4 3 

2 

1 

Abira R. 

.053 

.064 

.092 

.086 

.063 .135 .095 

.004 

.002 

.005 

.022 

.007 

.005 

.009 .008 

.005 

2 

Otsuchi Bay 

.081 

.094 

.120 

.113 

.087 .172 .123 

.013 

.003 

.003 

.014 

.002 

.014 

.006 .003 

— 

3 

Hirose R. 

.082 

.095 

.121 

.114 

.089 .169 .123 

.016 

.007 

.009 

.011 

.003 

.014 

.002 — 


4 

Kinu R. 

.069 

.080 

.104 

.097 

.074 .152 .110 

.014 

.009 

.011 

.015 

.008 

.011 

— 


5 

Kano R. 

.050 

.058 

.080 

.079 

.061 .103 .079 

.010 

.010 

.016 

.032 

.016 

— 



6 

Aga R. 

.083 

.094 

.123 

.116 

.090 .163 .115 

.016 

.005 

.007 

.005 

— 




7 

Iburibashi R. 

.080 

.090 

.113 

.104 

.081 .158 .108 

.025 

.018 

.016 

— 





8 

Yoshii R. 

.062 

.073 

.095 

.087 

.065 .152 .102 

.007 

.001 

— 






9 

Sufu R. 

.059 

.070 

.095 

.088 

.065 .146 .100 

.005 

— 







10 

Chikugo R. 

.031 

.040 

.062 

.054 

.035 .116 .076 

— 








11 

L. Biwa 

.035 

.028 

.052 

.045 

.045 .035 _ 









12 

Kumano R. 

.063 

.055 

.051 

.060 

.071 — 









13 

None R. 

.004 

.007 

.012 

.003 

— 









14 

Monobe R. 

.010 

.012 

.004 

— 










15 

Kure R. 

.022 

.025 

— 











16 

Shimanto R. 

.001 

— 












17 

Hitotsuse R. 

— 














the lower reach of each river were used for the 
comparison among populations. G S t value which 
indicates the extent of differentiation of allelic 
frequencies over 17 populations was estimated to 
be 0.544. 

D values which indicate the extent of differentia- 
tion between pairs of populations were shown in 
Table 5. D values between pairs of 17 populations 
ranged from 0.001 to 0.172. A dendrogram based 
on D values was further constructed to examine 
genetic relationships among the populations (Fig. 
3). Three major clustering groups were discrimi- 
nated on this dendrogram. Namely, the popula- 
tions from the Abira River to the Kano River in 
the Pacific coast and from the Aga River to the 
Chikugo River in the coast of the Sea of Japan 
(including the Yoshii River in the side of the 
Inland Sea), those from the None River to the 
Hitotsuse River in the Pacific coast, and those of 
the Lake Biwa and the Kumano River formed each 
of the three clustering groups. We gave names to 
these local groups; 1) Northern Group, 2) South- 
ern Group, and 3) Lake Biwa Group, respectively. 
From Table 5, D values between the Northern 
Group and the Southern Group ranged from 0.031 
to 0.123 (0.081 on average). Similarly, the values 


between the Northern Group and the Lake Biwa 
Group, and between the Southern Group and the 
Lake Biwa Group ranged from 0.076 to 0.172 
(0.125 on average), and from 0.028 to 0.070 (0.051 
on average), respectively. On the other hand, D 


— Kano R 
■ i— ChikugoR. 

- rSufu R. 

T Yoshii R. 
Abira R. 

— Iburibashi R. 
r A 9 aR - 

0 1 su chi B. 

rKinu R. 
'-HiroseR. 

L. Biwa 

Kumano R. 

| — KureR. 

|r None R. 




Monobe R. 
Shimanto R. 
Hitotsuse R. 


Northern 

Group 


Lake Biwa 
Group 


Southern 

Group 


005 0 

Genetic Distance 

Fig. 3. A dendrogram based on Nei’s genetic distance 
between pairs of 17 populations of T. hakonensis. 
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values between pairs of the populations within 
Northern, Southern, and Lake Biwa Groups 
ranged from 0.001 to 0.032 (0.010 on average), 
from 0.001 to 0.025 (0.010 on average), and 0.035, 
respectively. Thus, genetic distances between the 
local groups were considerably large compared 
with those within the groups. 

DISCUSSION 

Level of intrapopulational variability 

Intrapopulational variability at protein level in 
various species of animals has been examined and 
that of fishes is close to the standard levels in 
vertebrates [18]. In the populations of T. 
hakonensis, proportion of loci polymorphic and 
heterozygosity were estimated to be 0.22±0.10 
and 0.041+0.019, respectively (Table 3). These 
values are close to the mean values of Ostariophysi 
(0.16 and 0.045) and are within the range of fishes 
[18]. Therefore, it is regarded that intrapopula- 
tional variability of T. hakonensis is close to the 
standard levels in vertebrates. At the same time, 
this suggests that most of the populations in T. 
hakonensis has been maintained without extreme 
reduction of their size. 

Level of interpopulational differentiation 

Allelic frequencies were remarkably different 
among populations of three local groups in T. 
hakonensis (Table 3) and the extent of their 
differentiation was estimated by two kinds of 
indices. These data were further compared with 
those obtained in various animals. G ST value 
among the populations of T. hakonensis was 
estimated as 0.544. Nei [16] reported G ST value 
among populations in wild mice as 0.119, and that 
among populations in horseshoe crabs as 0.072, 
respectively. Thus, level of G ST value in T. 
hakonensis is several times as high as the levels 
hitherto reported. This reveals that most of 
genetic variation exists within populations in the 
animals hitherto examined, whereas a half of 
genetic variation exists among populations in T. 
hakonensis . 

Mean genetic distance between populations of 
different local groups in T. hakonensis was esti- 


mated as 0.051-0.125. Levels of genetic distance 
between local races or populations, and between 
subspecies of various vertebrates have been ex- 
amined hitherto. Nei [16] reported D values 
between local races and between subspecies in wild 
mice as 0.010-0.024 and 0.194, respectively. In 
fish, D values between populations were reported 
as 0.001-0.006 for salmonid fish [19] and as 
0.0007-0.0040 for red sea bream [20], and those 
between subspecies were reported as 0. 11-0.19 for 
cyprinid fish, the genus Campostoma [21]. Thus, 
level of D values between populations of local 
groups in T. hakonensis is much higher than the 
levels between local races or populations hitherto 
reported, and is close to those between subspecies. 
These results are consistent with those of 
mitochondrial DNA, since mitochondrial genomes 
had been also highly differentiated among three 
populations of each local group [22]. Level of D 
values between populations within each local 
group of T. hakonensis is a little higher than or 
close to the levels hitherto reported. 

In the other species of Tribolodon, T. brandti 
and T. ezoe y Sakai and Hamada [23] who ex- 
amined populations from the waters of Hokkaido 
pointed out that allelic frequencies at a few loci 
were largely different from populations of 
Fukushima Pref. we had reported previously [7]. 
Hanzawa [24] estimated the genetic distance be- 
tween these populations of two species based on 
allelic frequencies at 20 loci as 0.104 and 0.113, 
respectively. Thus, the high degree of interpopula- 
tionab differentiation is also documented in the 
other two species of Tribolodon. 

Geographical isolation of local groups 

The extent of genetic differentiation among the 
local groups of T. hakonensis was remarkably 
high. Fishes of each local group of T. hakonensis 
are very similar in morphology. Nakamura [25] 
who examined morphological characters of T. 
hakonensis could not clearly show the regional 
differentiation, although he found slight differ- 
ences among specimens collected from different 
regions in Japan. The present study on allozyme 
markers clearly showed the high degree of regional 
differentiation in the populations of T. hakonensis. 

We considered two posssibilities which could 
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explain the formation of these local groups. First is 
that any of populations of the local groups have 
undergone introgression and have accumulated 
genes derived from other related species. But, this 
is unlikely, because the local groups of T. 
hakonensis have been differentiated enough from 
the other species of Tribolodon at most of loci. 
Allelic differentiation between each local group of 
T. hakonensis and the other species of Tribolodon 
was observed at m-Aat, s-Idh-1 , 2, Ldh-1, 3 , Pgm, 
Sp-2 } 3 [6, 7, 24]. On the other hand, allelic 
differentiation among the local groups of T. 
hakonensis independently occurred at the other 5 
loci, Adh, a-Gdh-2 , Gpi-2, s-Mdh-2 , 6-Pgd. There- 
fore, we cannot interpret that the other related 
species had taken part in the interpopulational 
differentiation of T. hakonensis. 

The second possibility, which is the most prob- 


able, is that the populations of different local 
groups have been geographically isolated from 
each other since the time of the separation from 
their ancestral populations, and little gene flow has 
occurred among the populations of the groups. 
The dispersal of freshwater fishes is restricted by 
mountains, sea, etc., and regional differentiation 
occurs in their populations. Sakaizumi et al. [26] 
surveyed allozymic variation in Japanese wild 
populations of Oryzias latipes , and suggested that 
this case could be applied to interpopulational 
differentiation of this species because the bound- 
ary of the regions agreed well with watersheds or 
coastlines. 

Distribution of the three local groups of T. 
hakonensis was shown in Figure 4. Their distribu- 
tion shows distinct region specificity, and outline of 
their boundary agrees with watersheds or coast- 



Fig. 4. Distribution of three local groups of T. hakonensis. T. hakonensis has not been 
originally distributed in the area designated by g§. 
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lines. In the Southern Group of T. hakonensis, 
northeastern part of the boundary agrees with the 
Shikoku mountains, because T. hakonensis has not 
been originally distributed in the side of the Inland 
Sea of Shikoku ([27], Mizuno, N., personal com- 
munication). The Sikoku mountains may have 
restricted the dispersal of the Southern Group 
from the side of the Pacific Ocean to that of the 
Inland Sea. The restriction of dispersal by the 
Shikoku mountains is also suggested from the fact 
that faunae of freshwater fishes are largely dif- 
ferent between the side of the Inland Sea and that 
of the Pacific Ocean in Shikoku [27, 28]. Consider- 
ing that the Yoshii River population in the side of 
the Inland Sea of Honshu belongs to the Northern 
Group, the boundary of the Northern Group may 
agree with the coastline of the Inland Sea. In this 
regard, we speculate that T. hakonensis of the 
Northern Group has extended their distribution 
from the side of the Sea of Japan to that of the 
Inland Sea through plains around the Chugoku 
mountains. The reasons why the Northern Group 
had not extended their distribution to Shikoku 
through the Inland Sea is now obscure. But, we 
speculate one of the reasons is that T. hakonensis 
hardly migrates out of bays or coastal regions, as 
mentioned in detail later. Further investigation in 
the side of the Inland Sea of Honshu and Shikoku 
will give any suggestion. 

Lake Biwa group included two populations, but 
gene constitution of these populations was some- 
what different (Table 3, Fig. 3). Lake Biwa 
population which showed the unique gene con- 
stitution (Table 3) also suggests the geographical 
isolation. Lake Biwa had been originated in the 
Pliocene epoch, and has been maintained as the 
independent waters surrounded by watersheds 
[29]. Many endemic species, subspecies and local 
races are distributed in this lake, and it is con- 
sidered that they have been isolated within this 
lake and have been differentiated here. It is likely 
that T. hakonensis has been also isolated within 
the lake for a considerably long period of time and 
has been remarkably differentiated from the other 
populations. 

As stated above, the restriction of dispersal by 
watersheds may have affected the formation of the 
local groups. However, this cannot explain the 


reason why little gene flow is observed among the 
local groups although T. hakonensis can migrate 
through the sea regions. In this regard, we 
consider the following reasons. Though this spe- 
cies can migrate to the sea regions, it spawns in the 
middle reaches of rivers. This species is firmly 
dependent on the fresh waters, and may stay 
within bays or coastal regions without migrating to 
the open sea. Therefore, gene flow hardly occurs 
among populations of different local groups 
through the sea regions. This speculation is 
consistent with the observation that fish of this 
genus was distributed within the area of five meters 
in depth in Tokyo Bay [30]. Furthermore, this 
speculation may be supported by the observation 
that gene constitution of the amphidromous type 
of T. hakonensis was very similar to the landlocked 
type in the same waters (see Results). This 
observation suggests that the amphidromous type 
might have returned to its spawning ground in the 
middle reaches of the river without migrating to 
the other waters and had been mated with the 
landlocked type before these types were isolated 
by dums. 

It is likely that the populations of different local 
groups have been geographically isolated from 
each other for a long period of time. Divergence 
time between three local groups based on Nei’s 
genetic distance was estimated as approximately 
200,000-600,000 years, and agreed with the di- 
vergence time based on the net nucleotide differ- 
ences between mitochondrial DNA of three 
populations of each local group [22]. This time 
corresponds to the middle of diluvial epoch when 
the Japanese Islands had been strongly affected by 
ups and downs of the sea level with vicissitudes of 
glaciers and had connected with each other and the 
Asiatic Continent. It is likely that T. hakonensis 
had extended their distribution and that the local 
groups had been formed during this epoch. As for 
the dispersal and the differentiation of the local 
groups, we speculate as follows. Considering a fact 
that all of other related species of Tribolodon are 
distributed in the northern Japan and north of 
Japan at present, the ancestral population of T. 
hakonensis had been possibly distributed in or near 
these areas. First, the population of the Southern 
Group had extended their distribution to the 


460 


N. Hanzawa, N. Taniguchi and K. Numachi 


southern Japan and had begun to be genetically 
differentiated from the ancestral population. 
Secondly, a part of the population of the Southern 
Group had been isolated and the Lake Biwa group 
had been differentiated. On the other hand, the 
Northern Group had arisen from the ancestral 
population which had stayed in the north, and had 
extented its distribution rapidly. Further study as 
for specimens sampled at other waters will make 
sure this hypothesis. 
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